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Phenomenology of neutrino oscillations in vacuum and in cosmological plasma is 
considered. Neutrino oscillations in vacuum are usually described in plane wave 
approximation. In this formalism there is an ambiguity if one should assume 5p = 
and correspondingly SE 7^ or vice versa SE = and 5p ^ 0, or some other 
condition. We will use the standard set of quantum field theory rules and show 
that the theoretical description is unambiguous and that the state of oscillating 
neutrinos is determined by the conditions of their registration. The wave packet 
formalism is a natural result of such an approach. 

Oscillating neutrinos in cosmological plasma cannot be described in terms of wave 
function because of a fast loss of coherence due to elastic or inelastic scattering, 
so that one should use the density matrix formalism. Kinetic equations for the 
density matrix of oscillating neutrinos are derived. Numerical and semi-analytical 
solutions of the equations are examined. In particular, a possibility of amplification 
of cosmological lepton asymmetry in the sector of active neutrinos mixed to sterile 
ones is critically discussed. 

1 Neutrino oscillations in vacuum. Basic concepts. 

As is very well known, neutrinos possibly oscillate if (ancLhecause) their mass 
eigenstates are diffeijent from their interaction eigenstatestl □ (for a review and 
more references seeB). In other words, the mass matrix of different neutrino 
species is not diagonal in the basis of neutrino flavors: [i>e, f/i, i^t]- The latter 
is determined by the interaction with charged Icptons, so that a beam of e.g. 
electrons would produce v,, which is a mixture of several different mass states. 
An important condition is that the masses are different, otherwise oscillations 
would be unobservablc. Indeed, if all the masses are equal, the mass matrix 
would be proportional to the unit matrix which is diagonal in any basis. 

Of course not only neutrinos are capable to oscillate. All particles that 
are produced in the same reactions will do that, but usually the oscillation 
frequency, Wosc ^ 5m? I2E is so huge and correspondingly the oscillation length 

lose = 2p/5n? (1) 

is so small, that the effect is very difficult to observe. Here E and p are 



respectively the energy and momentum of the particles in consideration and 
Srri^ = nil ^ '^2- The expression is written in rclativistic limit. Only for Ja- 
mesons and hopefully for neutrinos the mass difference is so small that lose is, 
or may be, macroscopically large. 

The neutrino Lagrangian can be written as follows: 

jC^ =i9^i' + PMi' + P^i' + i?\fl (2) 

where the vector-column u = [vg, v^, Vt^" is the operator of neutrino field in 
interaction basis, I = [e, /x, t]'^ is the vector of charged lepton operators; the 
last two terms describe respectively neutral and charge current interactions 
(with Z and W bosons). The upper index "T" means transposition. The 
matrix Ai is the mass matrix and by assumption it is non-diagonal in the 
interaction basis. This is not necessary but quite natural because normally 
masses know nothing about interactions. 

Transformation between mass and interaction eigenstates is realized by an 
orthogonal matrix U with the entries that are parameters of the theory to be 
determined by experiment. In the simplest case of only two mixed particles 
the matrix U has the form 

U=( ^^'^M (3) 

y —smU cos J ^ ' 

If for example the only noticeable mixing is between electronic and muonic 
neutrinos, then the flavor eigenstates are related to the mass eigenstates z^i^2 
as: 

fe = Vx cos Q + sin d , 

f/j, = — I'l sin6' -|- 1^2 COS0 (4) 

Thus if electronic neutrinos are produced on a target by a beam of electrons, 
their propagating wave function would have the form 

V)^^(f,i) =cos6l|!^i)e''=i^+sin6l|i^2)e*'=^^ (5) 

where kx = ojt— kr and suh-Vg means that the initial state was pure electronic 
neutrino. Below (in this section only) we denote neutrino energy by (j to 
distinguish it from the energies of heavy particles that are denoted as E. We 
assume, as is normally done, plane wave representation of the wave function. 

If such a state hits a target, what is the probability of producing an elec- 
tron or a muon? This probability is determined by the fraction of fg and 



components in the wave function ip,y^ at space-time point x. The latter can be 
fomid by re-decomposition of 1/12 in terms of I'e^fj,'- 



^^^(f,i) =cos6le*Pi=" (cos6l|i/e) -sin6'|i'^)) +sin6le*P^^ (sin6l ji^e) + cose* ji^^)) 

One can easily find from that expression that the probability to register or 
i^fi are respectively: 

P^^ (r , t) - cos^ e + sin^ + 2 sin^ 9 cos^ 9 cos 5$ , (7) 



P^^ (f, - 2 sin^ 6* cos^ 6* (1 - cos 5$) (8) 

Here duj = lui ~ uj2, Sk ~ ki ~ k2, and 

(5$ (r, t) ^Scjt-Skr (9) 
The energy difference between the mass eigenstates is 

5ui = 5m^ + — 5k (10) 
dm? dk 

Using this expression we find for the phase difference 



= t + Sk \ -t-f] (11) 

The standard result of the neutrino oscillation theory are obtained if one 
assumes that: 1) Sk — 0, 2) k = ujv, and 3) t = r/v: 

<5$ = ^ (12) 
2k ^ ' 

Each of these assumptions is difficult to understand and, even more, some of 
them, in particular, 5k — Q may be explicitly incorrect (see below). Both the 
second and the third conditions are fulfilled for a classical motion of a point-like 
body, however their validity should be questioned for a quantum mechanical 
particle (for a wave). Despite all that the final result ( p^ is true and if there 
are some corrections, they can be trivially understood. 

Basic features of neutrino oscillatiDiia^were discussed in many papers. An 
incomplete list of references includes Qtj. One can find roare-jcitations and 
discussion in the above quoted papers and in the books Still some 

confusion and indications on possible controversies reappear from time to time 



in the literature, so it seems worthwhile to present a consistent derivation of 
eq. ( p^ ) from the first principles. A large part of this section is based on 
the discussions (and unpublished work) with A.Yu. Morozov, L.B. Okun, and 
M.G. Schepkin. 

Let us consider a localized source that produces oscillating neutrinos; we 
keep in mind, for example, a pion decaying through the channel tt ^ + v^. 
The wave function of the source ips i'r, t) can be Fourier decomposed in terms 
of plane waves: 



'Eot-ipoP I d^qC{q) exp -iq(f- Vot) = g^Sot-^Por^ y^^J (^3) 



where Vq — po/ Eq. It is the standard wave packet representation. The function 
C{p~po) is assumed to be sharply peaked around the central momentum 
Pq with dispersion Ap. The particle is, by construction, on-mass-shell, i.e. 
E'^ ~ p"^ + m? . This is of course also true for the central values Eq and 
Po- As the last expression shows, the particle behaves as a plane wave with 
frequency and wave vector given respectively by E^ and Pq and with the shape 
function (envelope) given by C{f— VqI) which is the Fourier transform of C{q). 
Evidently the envelope moves with the classical velocity Vq. Characteristic size 
of the wave packet is Ipack 1/ A.p. 

Let us consider the pion decay, n fi + v. One would naturally expect 
6u! ^ Sk ^ 5m? /E. If this is true the probability of oscillation would be 



Pose ^ cos 



x + b{x~Vt) 



lo 



(14) 



where lose is given by expression (|l|) and 6 is a numerical coefficient relating 
5p with lose- For simplicity the one-dimensional expression is presented. 

Thus to obtain the probability of neutrino registration one should aver- 
age the factor {x — Vt) over the size of the wave packet and for large pack- 
ets, if h is non-negligible, a considerable suppression of oscillations should 
be expected. The size of the neutrino wave packet from the pion decay at 
rest is macroscopically large, Ipack ~ ct.^ ~ 7.8 m, where c is the speed of 
light and Ttt = 2.6 • 10~* sec is the pion life-time. The oscillation length is 
lose = 0.4 m (i<^/MeV/((5m^/eV^), so lose could be smaller or comparable to 
Ipack and the effect of suppression of oscillations due to a finite size of the 
wave packet might be significant. It is indeed true but only for the decay of a 
moving pion, and this suppression is related to an uncertainty in the position 



of decay. To see that we have to abandon the naive approach described above 

and to work formally using the standard set of quantum mechanical rules. 

Let us assume that neutrinos are produced by a source with the wave 
function ips{x,t). This source produces neutrinos together with some other 
particles. We assume first the following experimental conditions: neutrinos 
are detected at space-time point Xv,tv, while the accompanying particles are 
not registered. The complete set of stationary states of these particles is given 
by the wave functions ijjn ^ exp(i_E„i). The amplitude of registration of prop- 
agating state of neutrino of type j (mass eigenstate) accompanying by other 
particles in the state i/'n is given by 



In principle one even does not need to know the concrete form of ipn: the 
only necessary property of these functions is the condition that they form a 
complete set: 



However in what follows we will use for simplicity the eigenfunctions of mo- 
mentum, '4>n ^ exp{ipr — iEt). 

For the subsequent calculations we need the following representation of 
the Green's function which is obtained by the sequence of integration: 



Here oj and p are respectively the fourth and space components of the 4- 
vector Pi. We have omitted spin matrices because the final result is essentially 
independent of that. The integration over dp was extended over the whole 
axis (from — oo to -(-oo) because the integrand is an even function oi p. This 
permit to calculate this integral by taking residues in the poles on mass shell: 
p = iVw^ — m? + ie. Both poles give the same contribution, so skipping 
unnecessary numerical coefficients, we finally obtain: 




(15) 




(16) 
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(17) 




-\-oo 



(18) 



As a source function ips we will take essentially expression ( p^ ) but assume 
that the source is a decaying particle with the decay width 7, that was born 
at the moment i = 0. It corresponds to multiplication of tps by 9{t) exp(— 7i). 
Taking all together we obtain the following expression for the amplitude: 

Jo J \r- rs\ J 

+00 

dwje''^'(*-*=)-''=^l^-^=l (19) 

-00 

Integrals over dfg and dp are taken over all infinite space. It is worthwhile 
to remind here that all the momenta are on mass shell, — + ra^ (we 
assumed that the source is a decaying pion) and lu"^ = k"^ + m^, where nij is 
the mass of j-th neutrino mass eigenstate. 

The integration over dts is trivial and gives the factor [E — En — LOj + 17/2) 
The integration over dfs can be easily done if the registration point is far from 
the source. In this case it is accurate enough to take l/|r — fs\ ~ l/r, while 
the same quantity in the exponent should be expanded up to the first order: 

\r^rs\~r~irs (20) 

where ^ = r/r is a unit vector directed from the center of the source taken 
at the initial moment i = to the detector at the point r. In this limit the 
integral over drg gives 5 (jj — pn — kj^, ensuring momentum conservation: 

p=Pn+kj=p^j (21) 

The vector of neutrino momentum is formally defined as 



fcj = ^ fcj W ^1 - (22) 



Ultimately we are left with the integral: 

^•+00 Auj^t—ik^r 



- / diVjC! (pn + kj - po) — -— - (23) 



where E^^j — y (p„ + fcj)^ + mj. This integral can be taken in the 'pole 
approximation' and to do that we need to expand the integrand near the 
energy conservation law, see below eq. (p^), as follows. The neutrino energy 
is presented as tOj = + Aujj. To avoid confusion one should distinguish 



between the deviation of neutrino energy from the central value given by the 
conservation law, Ao^j, from the difference of energies of different neutrino 
mass eigenstates, Slo — lo\ — 0^2- The neutrino momentum is expanded up to 



the first order in Awj : 



fc, - Juj] - m] « kf + Auj./Vj'^'^ (24) 



where vj-'^"' = k^^^ /uj^p is the velocity of j'-th neutrino. The pion energy is 
determined by the momentum conservation (^) and is given by 



^\l[Pn + kf + Afc,) + ml « 4°] + Aw, (25) 

where the pion velocity is vj'^'^ — (p„ + k'^^'^ ) / • The neutrino energy, wj*^^ , 
satisfying the conservation law is defined from the equation: 

Now the integral over ujj is reduced to 



4°] ~ - iof = (26) 



3 



r 

+00 giAcjj (t-r/V/) 

The last integral vanishes if i < vj'^\, while in the opposite case it can be 
taken as the residue in the pole and we finally obtain: 

4" - '-^^^ (^ - >f "') ("'°'' - "f^ - -2^f^^ (2") 

We obtained a neutrino wave packet moving with the velocity vj'^^ with a 
well defined front (given by the theta-function) and decaying with time in 
accordance with the decay law of the source. The similar wave packet, but 
moving with a slightly different velocity describes another oscillating state f^. 
It is evident from this expressions that phenomenon of coherent jpscillations 
takes place only if the packets overlap, as was noticed long agofflfl. 

The probability of registration of oscillating neutrinos at the space-time 
point (r, t) is given by the density matrix 

p,, = f dpn aP (f„ t,) (f,, t,) (29) 



The oscillating part of the probability is determined by the phase difference 
but now the quantities 5uj and 5k are unambiguously defined. To this end we 
will use the conservation laws (|2l|,|2^). They give: 



(5tj("^ = 5E^ and (5fc(°) = 5p^ (30) 
The variation of neutrino energy is given by 

5LU^V^''Hk + 5m^l2uj (31) 
while the variation of the pion energy can be found from expression (p5|): 

= yW^fc (32) 

From these equations follows 

(5cj = — ^ — ^ and 5k = ^ (33) 

One sees that generally both 5uj and 5k are non-vanishing. Only in the case 
of pion decay at rest, 5lo — Q but 5k is non-zero in any case. Substituting the 
obtained results into expression (O) for the phase difference we come to the 



standard expression (12) if = 0. This result shows a remarkable stability 
with respect to assumptions made in its derivation. However if the pion is 
moving, then the oscillation phase contains an extra term 

5<^ = ^ ' = — ; V — '-^—^ (34) 

2uj yM-^^W 2k 

This extra term would lead to a suppression of oscillation after averaging over 
time. This suppression is related to the motion of the source and reflects the 
uncertainty in the position of pion at the moment of decay. So this result can 
be understood in the frameworks of the standard naive approach. 

Similar expression can be derived for the case when both neutrino and 
muon from the decay n ^ fi + Vj are registered in the space^time points ?v, 
and Tfj.jtfj^ respectively. This case was considered in refs. EScJ. Here we will 
use the same approach as described above when the muon is not registered. 
The only difference is that in eq. ( p^ for the oscillation amplitude we have to 
substitute the Green's function of muon Gp(?Yi— r^, t^—tg) instead of '0n(^s, tg). 
The calculations are essentially the same and after some algebra the following 



expression for the oscillation amplitude is obtained: 



exp 



2{V^ + K) 

(fc(°)r, + fc(")r.-4°)i,-4")t.)] (35) 

where L = Vt — r. Each kinematic variable depends upon the neutrino state 
J, so they should contain sub-index j. The upper indices "0" mean that these 
momenta and energies are the central values of the corresponding wave packets, 
so that the classical relation — holds for them. Here the direction 

of momenta are defined as above along the vector indicated to the observation 
point. However, the kinematics in this case is different from the previous one 
and the change of the energy and momentum of each particle for reactions 
with different sorts of neutrinos are related through the equations: 

SE^ + 6E^ = and 5k^ + 5K = 0. (36) 

For the central values of momenta the following relations are evidently true, 
Ski, — 5Ei, — 5m? /2ki,, and the similar one for the muon without the last term 
proportional to the mass difference. Correspondingly the phase of oscillation 
is given by the expression 

(5$ = 5k^r^ - 6E^t^ + 6ki,ry - SE„ti, = 



,3.) 

The first two terms in the phase are proportional to the argument of C in 
eq. ( p5| ) and thus they give the contribution equal to the size of the the source, 
i.e. to the wave packet of the initial pion. If the latter is small (as usually the 
case) we obtain again the standard expression for the oscillation phase. 

2 Matter effects 

Despite extremely weak interactions of neutrinos, matter may have a signifi- 
cant infiuence on the oscillations if/because the mass difference between the 
propagating eigenstates is veoz. small. Description of neutrino oscillations in 
matter was first done in ref. Eil. Somewhat later a very important effect of 
resonance neutrino conversion was discoveredEa, when even with a very small 
vacuum mixing angle, mixing in medium could reach the maximal value. 
Hamiltonian of free neutrinos in the mass eigenstate basis has the form: 



0^(1.2) ^( El \ 
\ {) E2 J 



(38) 



where Ej = ^/p^ + m^. In the interaction basis Hm is rotated by the ma- 
trix (|: 



^(a,b) ^ ™(i.2) TT-i ^ f cos^ 9Ei+ sin^ 9 E2 g sin 9 cos 9 

un^ u I g sin 61 cos 61 sin^ 6* + cos^ 61 



(39) 



Here g — 6n{^ jlE and we returned to the more usual notation E for neutrino 
energy. 

The interaction Hamiltonian is diagonal in the interaction basis and if only 
first order effects in the Fermi coupling constant, G^, are taken into account, 
then the Hamiltonian can be expressed through refraction index, ria, of flavor 
a-neutrino in the medium (recall that the deviation of refraction index from 
unity is proportional to the forward scattering amplitude and thus contains 
Gi? in the first power): 

(^a,b) _( E(na-\) \ 

™* " E{n^~\) J 

where a small difference between Ei and E2 in front of small factors (n — 1) 
was neglected. 

Thus, up to a unit matrix, the total Hamiltonian in the interaction basis 
can be written as 

rT{a,b) _( I 5 sin 261/2 \ 
" 5 sin 20/2 ) 

where f = g cos 29 + ESn and Sn — Ua^ni,. This matrix is easy to diagonalize. 
Its eigenvalues are 



f±^P + g^sm^29 
Ai,2 - ^ (42) 

and the eigenstates in matter (up to normalization factor) are 

Refraction index may change with time, as happens in cosmology, or with 
space point, if neutrinos propagate in inhomogeneous medium, for example 
from the center of the Sun to its surface. If somewhere (or sometime) / vanishes 
then the resonance transition of one neutrino species to another is possiblecS. 
Indeed let us assume that i/g and are mixed with a small vacuum mixing 



angle 9 and that initially an electronic neutrino was produced in vacuum. So 
the initial propagating state would be mostly v,,: 

\vi)„, = \ue) + {l/2)t&n2e\u^) (44) 

After propagation in the media where the function / changes sign passing 
through zero, the propagating state would become mostly v^: 

\yl)f^n^\ye)~{2|ta.n2e)\u^,) (45) 

This effect of resonance conversion may play an important role in the resolution 
of the solar neutrino problem and in cosmology. 

3 Neutrino oscillations in cosmology 

3.1 A brief (and non- complete) review 

Neutrino oscillations in the primeval plasma is significantly different from e.g. 
solar neutrino oscillations in the following two important aspects. First, cos- 
mological plasma is almost charge symmetric. A relative excess of any particles 
over antiparticles is believed to be at the level 10~^ — 10"^", while in stellar 
material the asymmetry is of order unity. Neutrino oscillations in the early uni- 
verse may change the magnitude of asymmetry in the sector of active neutrinos. 
On the other hand, the asymmetry has a strong influence on the oscillations 
through the refraction index of the primeval plasma (see below). It leads to a 
strong non-linearity of the problem and makes calculations quite complicated. 
Second important point is that neutrino mean free path in the early universe 
is quite small at high temperatures and hence breaking of coherence becomes 
essential. Because of that one cannot use wave functions for description of os- 
cillations and should turn to the density matrix formalism. It leads to a great 
complexity of equations. Kinetic equations for density matrix with the arrnnnt 
of nnutrino scattering and annihilation were derived in the papers ErOoa. In 
ref. where the impact of neutrino oscillations on big bang nucleosynthesis 
(BBN) was first considered, only the second order effects, proportional to 
were taken into account, while the deviation of refraction index from unity 
was neglected. This approximation is valid for a sufficiently high (5m^. In a 
subsequent paper L3 implications for BBN of possible CP-violating effects in 
oscillaiinns. were discussed. Earlier works on neutrino oscillations also include 
refs.LjH'Lj. The^calculations of refraction index in cosmological plasma were 
performed in ref. cHI. The results of this work permitted to make_mQre accurate 
calculations of the role played by neutrino oscillations in BBNciJ lJ. 

It was noticed in ref.c3 that the oscillations between an active and sterile 
neutrinos could generate an exponential rise of lepton asymmetry in the sector 



of active neutrinos. The origin of this instabihty is the foUowing. Since lepton 
asymmetry comes with the oppo site sign to refraction indices of neutrinos and 
antineutrinos (see below section ^^), it may happen that the transformation 
of antineutrinos to their sterile partners would proceed faster than similar 
transformation of neutrinos, especially if resonance conditions are fulfilled. It 
would lead to an increase of the asymmetry and through the refraction index to 
more favorable conditions for its rise. However it was concludedc3 that the rise 
is not significant and the effect of the generated asymmetry on BBN is small. 
This conclusion was reconsidered in ref.Ej where it was argued that asymmetry 
generated by this mechanism could reach very large values, close to piity, and 
this effect, in accordance to the earlier paper of the same group E3, would 
have a significant influence on primordial abundances. This result attcacted 
a great attention and was confirmed in several subsequent publications E^ro. 
Moreover, some works shovpdrJapt only rising and large asymmetry but also a 
chaotic behavior of its signE3n^. 

Due to complexity of equations some simplifying approximations were 
made in their solutions. First -accurate numerical solution of (almost) ex- 
act equations were done in refs. eZI~0. The most essential approximation was 
that the lost of coherence was described by the term ^{peq — p) instead of the 
exact two-dimensional collision integral taken from the square of the scatter- 
ing amplitude over momenta. InjiacJ—this approximation was used in all the 
papers. The calculations of refs.EZHllI were performed for rather small mass 
difference, 6m? < 10^^ eV^. Neither chaoticity, nor a considerable rise of the 
asymmetry were found. For a larger mass difference a strong numerical insta- 
bility was observed. However this result is not in a contradiction with other 
papers because the latter found_.the above mentioned effects for much larger 
mass differences. In our workES we tried to extend the range of validity of 
direct numerical calculations to large Sm^. We were able to proceed only up 
to Sm? w 10~^. For higher values we did not find a way to avoid numeri- 
cal instability. To overcome the problem we analytically transformed kinetic 
equations to a simpler form that permitted a stable numerical integration. To 
this end an expansion in terms of a large parameter that corresponds to a high 
frequency of oscillations was used. This method is rather similar to the well 
known separation of fast and slow variables in differential equations. Accord- 
ing to our results the asymmetry may rise, but only by 5-6 orders of magnitude 
(from initial ^*^ ~^|^)| | n ^M''^ii"'i '^'^^'^ '^^ magnitude weaker than was obtained 
in other papers oE3e3'E3oE3i£j. However in agreemen±_with these papers no 
chaoticity was found. Our results were criticized in ref.Ej where it was argued 
that lepton asymmetry generated by oscillations must be chaotic. However the 
author misunderstood our calculations and criticized us for the things that we 



never did and, second, the approximation used in that paper E3, namely cal- 
culations in terms of one fixed value of neutrino momentum that was chosen 
to be equal to the thermally averaged one, is intrinsically inappropriate for 
the solution of the problem of chaoticity. In that approximation even neutrino 
oscillations in vacuum would result in chaotic lepton asymmetry. 

Thus, most probably a chaotic amplification of the lepton asymmetry does 
not take place, however the exact magnitude of the amplification remains un- 
certain. As discussed in ref. their recent calculations are exact and do 
not suffer from numerical instability. On the iither hand, we do not see any 
shortcomings of our semi-analytical approach l3, all approximations are well 
under control and numerical part of calculations is quite simple and stable. 
More work is necessary to resolve the contradictions both in magnitude and 
possible chaoticity. 



3.2 Refraction index 



In this section we derive the Schroedinger equation for neutrino wave function 
in the primeval plasma. We will start with the neutrino quantum operator 
h'aix) of fiavor a that satisfies the usual Heisenberg equation of motion: 



4 cos 9t, 



-Za(a;)OWz/a(x) = 0(46) 



where Ai is the neutrino mass matrix, W{x), Z{x) and e(x) are respec- 
tively the quantum operators of intermediate bosons and electrons, and 0„ — 
7a (1 ± 75). We assumed that the temperature of the plasma is in MeV range 
and thus only electrons are present in the plasma. 

Equations of motion for the field operators of W and Z bosons have the 
form 



G 



9 



2V2 



(47) 



i^a(x)Oi+Va(.T) + 



4 COS 9w 

(2 sin2 Ow - 1) e{x)Oi+'>e{x) + 2 sin^ dw e{x)Oi-^e{x) 



(48) 



where the differential operators G^ ^ are inverse Green's functions of W and 
Z. In momentum representation they can be written as 



It can be be shown that the term qaqp/rn^ gives contribution proportional to 
lepton masses and can be neglected. 

In the limit of small momenta, q ^ itiw,z, equation (E^ can be solved as 



9 



2^2 7 



1 



92 



'■w 



(e(a;)Oi+'z.e(x) 



(50) 



A similar expression with an evident substitution for the r.h.s. can be obtained 
for Za{x). These expressions should be substituted into eq. ( ^ ) to obtain 
equation that contains only field operators of leptons, i'a(a;) and e{x): 



2 



5") (e»0W^e(2;))l OWe(x) + 

Vb{x)0^^'^v^(x) + (2sin2 - l) e(x)0(,+)e(x) + 

2sin2 0H' e(:E)Oi-)e(:r))] Oi+V,(:r)} (51) 

Neutrino wave function in the medium is defined as 

vI/,(:r) = (A|z/,(a:)|A + z.W) (52) 

where A describes the state of the medium and v^^"* is a certain one-neutrino 
state, specified by quantum numbers fc, e.g. neutrino with momentum k. Equa- 
tion of motion for this wave function can be found from eq. ( |5l| ) after averaging 
over medium. The theory is quantized perturbatively in the standard way. We 
define the free neutrino operator vf^ that satisfies the equation of motion: 



(i^- M.)v^^\x) = 



(53) 



This operator is expanded as usually in terms of creation-annihilation opera- 
tors: 



J2[alu^{k)e^'^^+bl^v^{ky 



ikx 



(54) 



and one-particle state is defined as \v^''^) — a|,|vac). 

The equation of motion for the neutrino wave function '^a{x) can be ob- 
tained from expression ( ^ ) perturbatively by applying the operator (i^ — M.) 
and using eq. ( pT] ) with free neutrino operators in the r.h.s. After some 
algebra which mostly consists in using equations of motion for the free fermion 



operators and (anti)conimutation relations between the creation/annihilation 
operators, one would obtain the equation of the form: 



= + (55) 
where TLm is the free Hamiltonian; in the mass eigenstate basis it has the form 



Tio = diag + m?- . The matrix-potential Vg/ / describes interactions of 

neutrinos with media and is diagonal in the flavor basis. Up to the factor E (i.e. 
neutrino energy) it is essentially the refraction index of neutrino in the medium. 
The potential contains two terms. The first one comes from the averaging of the 
external current J lOal- Due to homogeneity and isotropy of the plasma only 
its time component is non- vanishing and proportional to the charge asymmetry 
(i.e. to the excess of particles over antiparticles) in the plasma. This term has 
different signs for neutrinos and antineutrinos. However the interactions of 
neutrinos with the medium is not always of the (current) x (current) form due 
to non-locality related to the exchange of W or Z bosons. If incoming and 
outgoing neutrinos interact in different space-time points, the interaction with 
the medium cannot be written as an interaction with the external current. The 
contribution of such terms is inversely proportional to ^ but formally it is 
of the first order in Gf- With these two types of contributions the diagonal 
matrix elements of the effective potential for the neutrino of flavor a has the 
form: 



K}/ = iCir/G^T^ + C^^3:I^ , (56) 



where E is the neutrino energy, T is the temperature of the plasma, Gp = 
1.166 • 10~^ GeV~^ is the Fermi coupling constant, a = 1/137 is the fine 
structure constant, and the signs "±" refer to anti-neutrinos and neutrinos 
respectively (this choice of sign describes the helicity state, negative for v 
and positive for D). According to ref . e3 the coefficients Cj are: Ci ~ 0.95, 
C| ~ 0.61 and 6*2'^ ~ 0.17. These values are true in the limit of thermal 
equilibrium, but otherwise these coefficients are some integrals from the dis- 
tribution functions over momenta. For oscillating neutrinos deviations from 
thermal equilibrium could be significant and in this case expression (^6|) should 
be modified. However it is technically rather difficult to take this effect into 



account in numerical calculations and the simplified version (56) is used. 



3.3 Loss of Coherence and Density Matrix 

Breaking of coherence appears in the second order in the Fermi coupling con- 
stant Gp, so that equations of motion for the operators of all leptonic fields 



(including electrons) should be solved up to the second order in Gp- Since the 
calculations are quite lengthy, we only sketch the derivation here. In the con- 
sidered approximation the lepton operators l{x), where I stands for neutrino 
or electron, in the r.h.s. of eq. ( |5l| ) should be expanded up to the first order in 
Gp- The corresponding expressions can be obtained from the formal solution 
of eq. ( IsTI ) up to first order in Gp- Their typical form is the following: 

l^la + Gi*{Y.h.s.o) (57) 

where the matrix (in neutrino space) Gi is the Green's function of the cor- 
responding lepton and r.h.s.o is the right hand side of eq. ( pl| ) taken in the 
lowest order in Gp, i.e with lepton operators taken in the zeroth order, I — Iq. 



Expression ( p7[ ) should be substituted back into eq. (51) and this defines the 
r.h.s. up to the second order in Gp in terms of free lepton operators Iq. Of 
course in the second approximation we neglect the non-local terms, ~ l/m^^. 

Now we can derive kinetic equation for the density matrix of neutrinos, 
Pj = i^V*, where over-hut indicates that p is a quantum operator. The C- 
valued density matrix is obtained from it by taking matrix element over the 
medium, p= (p). We should apply to it the differential operator {i^ — Ai) and 
use eq. (ply). The calculations of matrix elements of the free lepton operators 
Iq are straightforward and can be achieved by using the standard commutation 
relations. There is an important difference between equations for the density 
matrix and the wave function. The latter contains only terms proportional 
to the wave function, idf 'i' = TY^*, while equation for density matrix contains 
source term that does not vanish when p = 0. Neutrino production or destruc- 
tion is described by the imaginary part of the effective Hamiltonian. The latter 
is not hermitian because the system is not closed. By the optical theorem the 
imaginary part of the Hamiltonian is expressed through the cross-section of 
neutrino creation or annihilation. Such terms in kinetic equation for the den- 
sity matrix are similar to the "normal" kinetic equation for the distribution 
functions: 

II = I-oii (58) 

where the collision integral is the integral over the proper phase space from 
the following combination of the distribution functions 

F = -/l/2(l - /3)(1 - .A) + /4/3(l - /l)(l - /2) (59) 

A similar combination appears for the case of oscillating neutrinos but with a 
rather complicated matrix structure. For example there can be the contribu- 
tion of the form of the anti-commutator: 



- {pi,g{l - P3)5P2(1 - P4)g} + (inverse reaction) 



(60) 



and a few more of different structure that results in a very lengthy expres- 
sion. The matrix g describes neutrino interactions. It is diagonal in the flavor 
basis and has entries proportional to matrix elements squared of the relevant 
reactions. The Fermi-blocking factors (1 — p) appear when one takes matrix 
elements of neutrino operators over the medium in which neutrino occupation 
numbers may be non-zero. 

We will not uresent here the complete form of the equation. It can be found 
e.g. in the paperE3. Moreover in all the applications a "poor man" substitution 
has been done: all terms describing neutrino production or destruction were 
mimicked by 

-T{p-p,,) (61) 

where p^q is the equilibrium value of the density matrix, i.e. the unit matrix 
multiplied by the equilibrium distribution function 

= [exp(£;/T - C) + 1] (62) 

and the coefficient F is the reaction rate (see below). 



3.4 Oscillations and lepton asymmetry 

As we have already mentioned oscillations between active and sterile neutrinos 
may induce a significant lepton asymmetry in the sector of active neutrinos. 
Now we will consider this phenomenon in some more detail. Basic equations 
governing the evolution of the density matrix are: 

i(dt - Hpdp)paa = FQ{psa~ Pas)/'2.~ iTa{paa- feq) , (63) 

i{dt~ Hpdp)pss = -Foipsa - Pas)/'2 , (64) 

i{dt ~ Hpdp)pas WaPas + Foipss ~ Paa)/'^ - i'^lPas , (65) 

i{dt - Hpdp)psa = -WoPsa- Fq{Pss- Paa)/'2.-iTiPsa , (66) 

where a and s mean "active" and "sterile" respectively, Fq — Sm^ sm26/2E, 

Wo = Srr? cos 26/2E + V;^fj , H = ^^-KptotliM^ is the Hubble parameter, p is 

the neutrino momentum. The antineutrino density matrix satisfies the similar 
set of equations with the opposite sign of the antisymmetric term in VQ^-j and 
with a slightly different damping factor 7 (this difference is proportional to the 
lepton asymmetry in the primeval plasma). 

Equations ( |63[]6^ ) account exactly for the first order terms described by the 
refraction index, while the second order terms describing breaking of coherence 



are approximately modeled by the. damping coefScients Tj in accordance with 
eq. (|6ll). The latter are equal toO: 



OT^ 180C(3) 2 rp4 f(.'7\ 

To = 2Ti = ga GpT p . (67) 

In general the coefScient (7a (p) is a momentum-dependent function, but in the 
approximaiipn of neglecting [1 — /] factors in the collision integral it becomes 
a constants equal respectively to ~ 4 and gui_,,^i-r — 2.9l3. In the following 
we will use more accurate values found from the thermal average of the com- 
plete electro- weak rates (with factors [1 — £]-. included), which we calculated 
numerically from our Standard Model codella. This gives us g^^ ~ 3.56 and 
9y^4^^ - 2.5. 

It is convenient to introduce new variables 

X = rnQR{t) and y = pR{t) (68) 

where R{t) is the cosmological scale factor so that H = R/R and mo is an 
arbitrary mass (just normalization), we choose toq ~ 1 MeV. In the approx- 
imation that we will work, we assume that T — —HT, so that we can take 
R = 1/T. In terms of these variables the differential operator {dt — Hpdp) 
transforms to Hxdx- We will normalize the density matrix elements to the 
equilibrium function f^q: 

Paa = feq{y)[l + a{x,y)], Pss ^ feq{y)[l + S{x,y)] , (69) 

Pas = pI^^ feq{y)[h{x,y)+il{x,y)] , (70) 

and express the neutrino mass difference 6m? in eV^. 

As the next step we will take the sum and difference of eqs. (|63|)-(|6^) for 
v and 9. The corresponding equations have the following form: 

4 = Fl± > (71) 
a± = -Fl± - 27+a± - , (72) 

/i± = Ul±~VZl^~-f+h±~-f^h^ , (73) 

l± = ^{a±-s±)-Uh± + VZh^-j+l±-j_l^ , (74) 

where a± = (a ± a)/2 etc, and the prime means differentiation with respect 
to X. We have used W = U ± VZ, 7 = Ti/Hx, and 7± = (7 ± 7)/2, where 
7„ parameterizes the difference of interaction rates between neutrino and anti- 
neutrinos, which is proportional to the neutrino asymmetry. With the approx- 



imation ptot — 10.757r^T'*/30, the expressions for t/, y, and Z become: 

2 

U = 1.12 -10^ cos 26*^— + 26.24-, 

y x'^ 



where ?7o is the asymmetry of the other particle species: 

77o = 2iq^^ + 77^^ + 77^^ + 77e - 77,1/2 (for v^) , (78) 
€ = 277,^ + 77,^ + 77,^ - 77„/2 (for ly^) , (79) 

and r] for v^- is obtained from eq. ( |79| ) by the interchange /i ^ r. The asym- 
metry is normahzed in the same way as the neutrino asymmetry (the second 
term in (|7^)). Here we have imphcitly assumed that Va = v^. 

Up to this point our equations are essentially the same as those used 
by other groups. The equations look rather innocent and at first sight one 
does not anticipate any problem with their numerical solution. However the 
contribution from Z could be quite large with the increasing magnitude of 
the asymmetry. The exact value of the latter is determined by a delicate 
cancellation of the contributions from all energy spectrum of neutrinos. The 
function a_ under momentum integral is quickly oscillating and very good 
precision is necessary to calculate the integral with a desired accuracy. Even 
a small numerical error results in a large instability. To avoid this difficulty 
we analytically separated fast and slow variables in the problem and reduced 
this set of equations to a single differential equation for the asymmetry that 
can be easily numerically integrated. The corresponding algebra is somewhat 
complicated and we will not discuss it here. One can find the details in ref.lIS. 

We have found that asymmetry practically does not rise for a large mixing 
angles, sin 6* > 0.01. For smaller mixings some rise of the asymme±rv_is ob- 
served, though much weaker than that obtained in refs. oEjQ'ESEjEj'EI. For 
example for Sm^ — 1 eV^ wc found that asymmetry rises by approximately 
5 orders of magnitude reaching the value around 10^^ for sin 20 in the inter- 
val 10^^ — 10^'^. For Sm^ = 10^ eV'^ the symmetry could rise up to 0.01 for 
sin 26* = 3 • 10~^ — 3 • 10~^, and only for huge mass difference Sm^ = 10^ 
the asymmetry may reach unity in a rather narrow range of mixing angles, 
sin26' = 2 • 10~^ — 4 • 10~^. For sin 26* outside the indicated limits the asym- 
metry does not rise. For even larger mass difference, Sm'^ = 10^^ eV^, the 
asymmetry practically does not rise. 



(75) 
(76) 
(77) 



Resolution of the contradiction between different groups is very impor- 
tant for the derivation of the constraints on the parameters of the oscillations 
from big bang nucleosynthesis (BBN). There are several effects by which the 
oscillations may influence abundances of light elements: 

1. If sterile neutrinos are excited by the oscillations then the effective num- 
ber of neutrino species at nucleosynthesis would be larger than 3. This 
effect, as is well known, results in an increase of mass fraction of helium-4 
and deuterium. 

2. Oscillations may distort the spectrum of neutrinos and, in particular, 
of electronic neutrinos. The sign of the effect is different, depending 
on the form of spectral distortion. A deficit of electronic neutrinos at 
high energy results in a smaller mass fraction of helium-4, while a deficit 
of Vf. at low energy works in the opposite direction. A decrease of to- 
tal number/energy density of v,, would result in an earlier freezing of 
neutrino-to-proton ratio and in a larger fraction of helium-4. 

3. Oscillations may create an asymmetry between v^, and anti-i^g. The n/p 
ratio in this case would change as n/p ~ exp(/x/T), where \x is the chem- 
ical potential corresponding to the asymmetry. With the present day 
accuracy of the data the asymmetry in the sector of electronic neiitrinos 
could be at the level of a few per cent, i.e. much larger than the standard 
10~^°. Even if asymmetry is strongly amplified but still remains below 
0.01 its direct influence on BBN would be negligible. It may however 
have an impact on the nucleosynthesis in an indirect way. Namely, the 
rise of the asymmetry by several orders of magnitude could suppress neu- 
trino oscillations through refraction index so that new neutrino species 
corresponding to sterile neutrinos arc not efficiently excited. 

Thus one sees that the effects of oscillations may result both in a reduc- 
tion or in an increase the effective number of neutrino species. In particular, 
oscillations may open room for additional particles at BBN. This conclusion 
strongly depends upon the magnitude of lepton asymmetry generated by the 
oscillations. Thus a resolution of the controversies between different theoretical 
calculations would be very important. 
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